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Abstract 

A direct hnk between the topological complexity of ferromagnetic media 
and their dynamics has recently been established through the construction of 
unambiguous conservation laws as moments of a topological vorticity. In the 
present paper we carry out this program under completely realistic conditions, 
with due account of the long-range magnetostatic field and related boundary 
effects. In particular, we derive unambiguous expressions for the linear and angular 
momentum in a ferromagnetic film which are then used to study the dynamics of 
magnetic bubbles under the influence of an applied magnetic-fleld gradient. The 
semi-empirical golden rule of bubble dynamics is verifled in its gross features but 
not in its flner details. A byproduct of our analysis is a set of virial theorems 
generalizing Derrick's scaling relation as well as a detailed recalculation of the 
fundamental magnetic bubble. 
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I. Introduction 



Magnetic bubbles have been known to exhibit some distinct dynamical 
features due to their nontrivial topological structure. The inherent link between 
topology and dynamics was already apparent in the early work of Thiele [1] as well 
as in many investigations that followed [2, 3]. The essence of the early work is best 
summarized by the experimentally observed skew deflection of magnetic bubbles 
under the influence of an applied magnetic-field gradient. The so-called golden 
rule of bubble dynamics relates the deflection angle 6 to the winding number Q by 

2 

|^sin5 = Q, (1.1) 

where g is the strength of the applied field gradient, r is the bubble radius and V 
its speed. Relation (1.1) is remarkable in two respects. First, it suggests that only 
topologically trivial {Q = 0) bubbles move in the direction of the gradient (6 = 0), 
even though such a behavior would naively be expected for all magnetic bubbles; 
in fact, bubbles with a nonvanishing winding number {Q = ±1, ±2, . . .) tend to 
be defiected in a direction nearly perpendicular {S ~ 90°) to the applied gradient. 
Second, Eq. (1.1) implies some sort of a topological quantization in that it relates 
the integer-valued winding number to experimentally measured quantities that 
can, in principle, assume any values. 

Although golden rule (1.1) has been employed with considerable success in 
the analysis of actual experiments, especially for hard (|Q| » 1) bubbles, it is 
still a semi-empirical relation whose precise meaning and domain of validity must 
be specified. For instance, the meaning of the various quantities entering Eq. 
(1.1) needs to be explained because magnetic bubbles are extended structures 
rather than point-like particles. Furthermore the usual derivations of (1.1) are 
based on the assumption that the bubble reaches a steady state, in the presence 
of dissipation, in which the deflection angle 5, the radius r and the speed V 
approach constant values. But such an assumption was never justifled and is 
actually incorrect. In practice, experiments are analyzed by applying Eq. (1.1) 
with average values for the deflection angle and the speed and by assuming that 
the radius does not change signiflcantly during the application of the gradient. 

In some recent work on this subject [4] the link between topology and 
dynamics was made explicit through the construction of unambiguous conservation 
laws as moments of a suitable topological vorticity. The important qualitative 
features of bubble dynamics became then apparent. Thus, in the absence 
of external magnetic-fleld gradients or other perturbations, bubbles with a 
nonvanishing winding number cannot move freely but are always spontaneously 
pinned. On the other hand, in the absence of dissipation, a bubble would be 
deflected at a right angle (5 = 90°) with respect to an applied magnetic-fleld 
gradient, with a drift velocity that can be calculated analytically in some important 
special cases [4] and is generally consistent with Eq. (1.1). The emerging picture 
is thus analogous to the Hall motion of an electron as well as to the Magnus effect 
of fluid dynamics. These analogies further suggest that the deflection angle should 
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deviate from 90° in the presence of dissipation. However an exact calculation of 
the deflection angle, i.e., a derivation of the golden rule, is no longer possible on 
the basis of conservation laws alone. 

Therefore the semi-quantitative picture derived from the conservation laws 
must be supplemented by some results from a numerical solution of the underlying 
Landau-Lifshitz equation. Such a solution is not straightforward under completely 
realistic conditions; calculation of the long-range magnetostatic field is always a 
problem and the finite thickness of actual magnetic films forces one to work with 
a three-dimensional (3D) grid, even though the essential topological structure of 
magnetic bubbles is two-dimensional (2D). Hence our numerical efforts have thus 
far been restricted to strictly 2D models of increasing complexity [5, 6]. 

It is the purpose of the present paper to provide a precise formulation within 
the quasi-2D geometry of realistic ferromagnetic films, taking into account the 
effects from the film boundaries and the magnetostatic field. In Section II 
we review some basic facts about the Landau-Lifshitz equation and introduce 
convenient (rationalized) physical units. In Section III we discuss the two 
ingredients that are important to establish a direct link between the topological 
complexity of magnetic structures and their dynamics, the gyrovector and the 
stress tensor. The derivation of unambiguous conservation laws is then carried 
out in Section IV in the presence of film boundaries. A byproduct of our study 
of conservation laws is a set of virial theorems that generalize the well-known 
scaling relation of Derrick [7] , an issue taken up in Section V. A detailed numerical 
calculation of the fundamental ((5 = 1) bubble is presented in Section VI which 
is consistent with the virial theorems. The issue of skew defiection in an applied 
magnetic-field gradient is studied in Section VII where we find that the semi- 
empirical golden rule is verified in its gross features but not in its details. Our 
conclusions are summarized in Section VIII together with some suggestions for 
further work. 



II. The Landau-Lifshitz equation 

A ferromagnetic medium is described in terms of the density of magnetic 
moment or magnetization M which is due primarily to the electron spins but 
may include contributions also from the orbital motion. In general, the vector 
M = (Ml, M2, M3) is a function of position and time except that its magnitude is 
nearly constant for a wide temperature range sufficiently below the Curie point. 
Thus we write 

M = M(x, t), = + M| + M| = M^, (2.1) 

where x = (xi, X2, xs) is the position vector, t is the time variable, and the constant 
Mo is the saturation magnetization. 
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Static as well as dynamical properties of the magnetization are governed by 
the Landau-Lifshitz equation 



+ 7(MxF) = — Mx— - , (2.2) 



dt ' M V dt 

which describes precession around an effective field F with the constant 7 given 

by 

7 = (2.3) 

zrrieC 

where ~ 2 is the gyromagnetic ratio, e the electron charge, me the 
electron mass, and c the velocity of light. Equation (2.2) also includes 
a phenomenological (Landau-Gilbert) dissipative term where the dissipation 
constant A is dimensionless. This choice of dissipation preserves the magnitude of 
magnetization. 

The effective field F may be written as 

F = Fe + F, + H6 + H. (2.4) 

Here Fg is the exchange field 

2A 

F. = ji^AM, (2.5) 

where A is the exchange stiffness constant and A the Laplace operator. F^ is the 
anisotropy field 

2K 

Fa = -j^{M„M„0), (2.6) 

where K is a positive constant leading to an easy axis in the third direction. In 
ferromagnetic films made out of bubble materials the easy axis is perpendicular 
to the film surface [2] . H5 is a uniform bias field, 

H6 = (0,0,if6), Hb = const, (2.7) 

applied along the easy axis. Finally H is the magnetic field produced by the 
magnetization itself and thus satisfies the magnetostatic equations 

VxH = 0, V-B = 0; B = H + 47rM, (2.8) 

where B is the corresponding magnetic induction. The use of the magnetostatic 
instead of the complete Maxwell equations is justified by the fact that time 
variations of magnetic structures of practical interest are slow. 

Numerical values of the various constants introduced above may be found in 
Ref. [2] for a number of ferromagnetic materials. However using all these constants 
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within a theoretical development clouds the underlying simplicity of the Landau- 
Lifshitz equation. Hence we introduce rationalized physical units as follows. First, 
we work with the normalized magnetization 

m = M/Mo, = 1. (2.9) 

Second, we measure distance, time and magnetic field (induction) in units of 



A 1 

-77^, z 7-r and AttMq, (2.10) 

27rM2' 47r7Mo ^ ^ 

respectively. We further define the dimensionless anisotropy constant 

which is usually referred to as the quality factor. Finally, we introduce new symbols 
for dimensionless magnetic fields, such as h = H/AnMo, but maintain the same 
symbols x and t for the rationalized space and time variables. 
The Landau-Lifshitz equation is then written as 

rii + (m X f) = A(m x m), = 1, (2.12) 

where the dot denotes diff'erentiation with respect to time, a convention that will 
be adopted from now on. The effective field f is given by 

f = Am - /t(mi,m2,0) + hfc + h, (2.13) 

where hf, = (0, 0, hb) is the bias field and h satisfies the magnetostatic equations 

Vxh = 0, V-b = 0; b = h + m. (2.14) 

The only free parameters are now the quality factor k, the dissipation constant A, 
and the bias field hi,. One should add that Eq. (2.12) differs from the Landau- 
Lifshitz equation studied in our earlier work [4] by the overall replacement m 
— m that originates in a sign difference between magnetization and spin introduced 
by the negative electron charge. The resulting sign differences at various stages of 
the calculation will be incorporated without further notice. 

The remainder of this section is devoted to a discussion of the hamiltonian 
structure associated with the Landau-Lifshitz equation at vanishing dissipation 
(A = 0). Then we write 

m + (m x f) = 0, = 1, (2.15) 

where the effective field f is still given by Eq. (2.13) and may be expressed entirely 
in terms of the magnetization once the magnetostatic field h is determined by 
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solving the linear system (2.14). Now a conserved energy functional W = W{m) 
exists such that the effective field is obtained through the general relation 



f = (2.16) 
dm 

where the symbol d denotes the usual functional derivative. Equation (2.16) 
together with Eq. (2.15) imply that the functional W is indeed conserved, and 
that (2.15) is the Hamilton equation associated with the hamiltonian W endowed 
with the Poisson-bracket relations 

{mi(x), mj(x')} = £yfcTOfc(x)5(x - x'), (2.17) 

which are reminiscent of the spin commutation relations. Here eijk is the 3D 
antisymmetric tensor and the usual summation convention for repeated indices is 
invoked. 

In order to display the explicit form of the energy functional we write 

W = We + Wa + Wb + Wm, (2.18) 

where the four terms correspond to the exchange, anisotropy, bias, and 
mangetostatic field. The more or less obvious choice of the exchange energy. 



J WedV, We^^idim-diUi), (2.19) 



where We is the corresponding energy density, requires some qualification in the 
presence of boundaries. Thus we consider the functional variation 

SWe = J {diSia ■ dim)dV = j{5m ■ diia)dSi - j {5m- Am)(iF, (2.20) 

where the surface-element vector dS = {dSi,dS2,dSs) is perpendicular to the 
boundaries of the ferromagnetic medium. Equation (2.20) would yield the desired 
relation 

Fe = = Am, (2.21) 

dm 

if the surface integral were absent; that is, if the gradient of the magnetization 
along the normal to the surface vanished. We write symbolically 

^ = 0, (...) 

which will be viewed as a boundary condition to be imposed at the free boundaries 
of the medium, in addition to the familiar boundary conditions of magnetostatics. 
This "unpinned" boundary condition was previously employed in the study of 
ferromagnetic films [2] and will play an important role in the following. 
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On the other hand, the usual bulk expressions for the anisotropy and bias 
(Zeeman) energies, 



W^a = y WadV, Wa = ^(^1 + ml) 

(2.23) 

Wb= WbdV, Wb = hb{l - ma). 



are free of boundary ambiguities and obviously yield the corresponding 
contributions to the effective field through the general relation (2.16). Note that 
in the Zeeman energy we have subtracted the (trivial) contribution from the state 

mo = (0,0,1), (2.24) 

which describes a fully saturated ferromagnet and will thus be referred to as the 
ground state. Configuration (2.24) is the simplest example of a static solution of 
the Landau-Lifshitz equation. 

Boundary effects are also important in the definition of the magnctostatic 
energy and will be described here in some detail. We begin with the reasonable 
ansatz 

Wm=l f h^dV, (2.25) 



where it is understood that the integral extends over all volume, inside and outside 
the material, and that the field h is expressed in terms of the magnetization 
through Eqs. (2.14). However, in order to justify that (2.25) is an appropriate 
choice of the mangetostatic energy within the context of the Landau-Lifshitz 
equation, one must show that 

h=-^. (2.26) 
dm 

Such a demonstration is not completely straightforward because of the implicit 
dependence of h on the magnetization. 

To make this dependence explicit we introduce a scalar potential from 

h=-VV', A^ = (V-m), (2.27) 

and solve the Poisson equation to obtain 



m(x') . dS' f ( V.m)(x') ^^^, 



X — X 



(2.28) 



where the surface integral extends over the boundaries of the ferromagnetic 
medium, if any, and the volume integral over the bulk of the medium. Applying a 
careful partial integration yields the equivalent relation 

1 /" (x — x') ■ m(x') , 
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whose advantage is that it contains no derivatives of the magnetization and is vahd 
irrespectively of the presence of boundaries. 

As an elementary illustration we consider a ferromagnetic film of thickness 
d (see Fig. 1) and assume that the magnetization is equal to its ground-state 
value (2.24) inside the film (region I) and vanishes outside (regions II and III). An 
explicit calculation of the integral in Eq. (2.29) then yields 

xs > d/2 

-d/2 <X3< d/2 (2.30) 
X3 < -d/2. 

Therefore the magnetic field is given by h = — VV' = —mo inside the film and 
vanishes outside. The magnetic induction b = h + mo vanishes everywhere. 

We now return to the magnetostatic energy (2.25) and replace by — h-(V'0). 
An application of the divergence theorem and Eqs. (2.14) then gives 




ijjim-dS) - / •0(V-m)dV 



(2.31) 



where we have also used the fact that t/^ is continuous across the boundary and 
that the difference between the normal components of the magnetic field on the 
two sides of the boundary is equal to the normal component of the magnetization. 
A further partial integration transforms (2.31) into 



m 



= ~J{h-m)dV, (2.32) 



which shares with Eq. (2.29) the property that it is valid whether or not 
boundaries are present. Hence, using this form of the magnetostatic energy and a 
magnetic field calculated from Eq. (2.29), the basic relation (2.26) is established 
by straightforward manipulations. 

We complete the discussion of the canonical structure noting that the Landau- 
Lifshitz equation is actually a constrained hamiltonian system. Nevertheless one 
may resolve the constraint = 1 explicitly using, for example, the spherical 
parametrization 

mi=sin0cos$, m2 = sin0sin$, 777.3 = cos©. (2.33) 

The energy functional is then parametrized in terms of the two independent fields 
© and $ and the general form of the Landau-Lifshitz equation reads 

SW ■ SW 

sin©© = - — , sin©$ = — , (2.34) 

which suggests that the pair of fields 

n = cos© and $ (2.35) 
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is a canonical pair: 

■ 5W ^ SW 

However most of the special dynamical features of the ferromagnetic continuum 
emerge exactly when the definition of the canonical variables (2.35) encounters 
ambiguities due to a possibly nontrivial topological structure of the magnetization. 

Finally we return briefly to the issue of dissipation and rewrite Eq. (2.12) in 
the equivalent form 



m + (m X G) = 0, = 1 



1 , A 



(2.37) 



G = Aif+A2(mxf), Ai = — — ^, A: 



1 + A2' " 1 + A2 



We then examine the rate at which the energy changes in the presence of 
dissipation: 

W = J ■ihjdV = -Aa J [f^ - (m • ff]dV. (2.38) 

Because m is a unit vector the integrand in the last step of Eq. (2.38) is positive 
definite and the energy decreases when the dissipation constant A is positive. 



III. Gyrovector and the stress tensor 

The key quantity for the description of both topological and dynamical 
properties of the magnetization is the gyrovector or vorticity = (71, 72, 73) whose 
cartesian components are given by 

7i = -^^ijkidjm X dkin) ■ m. (3.1) 

The former terminology was introduced in the early work [1] but the latter seems 
more appropriate in view of the significant formal analogy of the vector with 
ordinary vorticity in fiuid dynamics. Nevertheless one should stress that is not 
related to actual rotational motion in the ferromagnetic continuum but rather 
to the topological complexity of the magnetization. For the moment, we are 
concerned with instantaneous properties of the unit vector field m = m(x) at 
some instant t that is not displayed explicitly. Questions of dynamics will be 
addressed later in this section. 

An immediate consequence of the definition (3.1) and the constraint = 1 
is that the vorticity field is solenoidal, 

V - =0, (3.2) 
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and hence the corresponding vortex hnes cannot terminate but at the boundaries 
of the ferromagnetic medium. The precise nature of vortex hnes is revealed by 
expressing the vorticity in terms of the canonical variables (2.35), 

= Vn X V$, (3.3) 

a relation that suggests an analogy of 11 and $ with the Clebsch potentials of fluid 
dynamics [8] . It also establishes that vortex lines are defined as the intersections of 
the two surfaces n(x) = ci and $(x) = C2 where ci and C2 are arbitrary constants. 
In other words, vortex lines are the curves along which the magnetization vector 
m remains constant. 

Such a simple definition of vortex lines allows a transparent topological 
classification of the possible distributions of magnetization. We shall consider 
the physically interesting class of configurations m = m(x) that are differentiable 
functions of position and approach the ground state of the ferroniagnct at spatial 
infinity: 

m(x) ^ mo = (0,0,1). (3.4) 

X|^oo 

In the absence of boundaries the medium extends to infinity in all directions and 
vortex lines are closed curves. One may then define a degree of knottedness, or 
helicity, of tangled vortex lines by analogy with related work in fluid dynamics [9] 
and magnetohydrodynamics [10]. The current status of the topological aspects of 
the above subjects may be traced from Ref. [11]. In the present context, such 
a degree is more appropriately referred to as the Hopf index [12]. In view of 
the boundary condition (3.4) the 3D space is isomorphic to the sphere and a 
specific configuration m = m(x) establishes a map from to S'^, where S'^ is the 
2D sphere defined from the constraint = 1. Such a map is characterized by 
the integer-valued Hopf index defined as follows. Let m(x) = mi and m(x) = m2 
be any two vortex lines where mi and m2 are constant unit vectors. The linking 
number of these two curves is independent of the specific choice of the pair of 
vortex lines and is called the Hopf index of configuration m(x). 

In order to make a first contact with the dynamics we also quote an analytical 
definition of the Hopf index. The solenoidal vorticity is derived from a vector 
potential a, 

= V X a, (3.5) 

and the Hopf index is given by 

= ^ /(a • )dV. (3.6) 

Although the vector potential is unique only to within a gauge transformation, 
is gauge invariant and may be expressed entirely in terms of by 

^ = (4^ / ^-■^7.(x) j^^7.(x')^V^ciF'. (3.7) 
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The remarkable fact is that the above integral is always equal to an integer, and its 
explicit values coincide with those obtained through the linking-number definition 
given in the preceding paragraph [13, 14]. 

It would appear that a simpler (local) expression for the vector potential may 
be derived from Eq. (3.3) which suggests that 

a = nV$, (3.8) 

provided that 11 and $ are differentiable functions of position. However inserting 
(3.3) and (3.8) in (3.6) would then lead to a vanishing Hopf index. Putting it 

differently, the canonical variables 11 and $ cannot be both differentiable for field 
configurations with ^ 0, even though the magnetization is always assumed to 
be differentiable. Indeed explicit examples worked out in the literature [4, 15] 
demonstrate that the magnetization reaches the north as well as the south pole 
of the sphere = 1 along certain (vortex) lines where the angular variable $ 
becomes multivalued when N ^ 0. While these difficulties are largely irrelevant, 
because of the gauge-invariant definition (3.7), they already provide an important 
hint concerning dynamics. Note that the vector potential (3.8) coincides with 
the familiar expression for the momentum density associated with the Hamilton 
equations (2.36). We thus conclude that the difficulties discussed in connection 
with Eq. (3.8) render ambiguous also the usual linear momentum. 

We defer for the moment further discussion of dynamics and return to the 
question of topological classification in the presence of boundaries. Specifically we 
consider the film geometry of Fig. 1 where vortex lines need not be closed but 
may terminate at the boundaries of the film. Hence a definition of a Hopf index 
is no longer meaningful. Instead we consider the flux of vorticity 

g = i^/-dS (3.9) 



through any open surface 5" that is contained within the film but extends to infinity 
on all sides. The fiux is independent of the specific choice of a surface with the 
above properties, thanks to V- = and an elementary application of the divergence 
theorem. In particular, S may be a plane perpendicular to the third axis, 

<5 = ^ / lzdxidx2, --<Xi< -, (3.10) 

where the double integral is independent of 0:3. In fact, this integral coincides 
with the Pontryagin index or winding number [13] of the magnetization and is also 
integer-valued [Q = 0, ±1, ±2, . . .). Again, when Q 7^ 0, the canonical variables H 
and $ cannot be defined everywhere and the corresponding linear momentum is 
ambiguous. 

To be sure, the topological classification described above does not assume 
that the configuration m = m(x) solves the Landau-Lifshitz equation but merely 
that it obeys some general physical restrictions such as differentiability and Eq. 
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(3.4). Of course, this classification would become especially relevant if stationary 
solutions were found with a nontrivial topology. In this respect, we note that 
arguments of varying completeness have been presented in the literature for the 
existence of magnetic vortex rings with a nonvanishing Hopf index [4, 15] but a 
definite theoretical treatment and actual observation are still lacking. Nonetheless 
magnetic bubbles with a wide range of winding numbers have been observed in 
ferromagnetic films [2, 3]. Finally, if the magnetization itself is allowed to be 
nondifferentiable at isolated singular points, one is naturally led to a class of 
topological defects that are also characterized by a winding number of the form 
(3.9) except that the surface S is now closed around a singular point. Such defects 
have been observed in the bulk of the ferromagnetic continuum and are called 
Bloch points [2, 3]. 

We now organize the various hints concerning the connection between 
topology and dynamics by considering the time evolution of the vorticity (3.1). 
An elementary calculation based on the Landau-Lifshitz equation at vanishing 
dissipation, Eq. (2.15), leads to 

7i = -eijkdj{f- dkm.) = SijkdjTk, (3.11) 

where 

Tk = -(f • dkui) = f ^ ■ dkmj (3.12) 

is the "generalized force density" that appeared first in the work of Thiele [1]. We 
take this calculation one step farther using the formal argument 



5W \ 

— ■ dkm. dV = dkW = (3.13) 
dm / 

to conclude that may be written as a total divergence, 

Tk = deake, (3.14) 

where ak£ will be called the stress tensor. Equation (3.11) then reads 

ii = SijkdjdeCTke (3.15) 

and proves to be fundamental for our purposes [4]. 

To complete this line of reasoning we must also supply an explicit expression 
for the stress tensor. As a first step we insert in Eq. (3.12) the effective field f of 
Eq. (2.13): 

Tk = Tk+rk +Tk+T^, 

Tk = -(Am • dkm), Tk = K{midkmi + m25fc"^2), (3.16) 
= -hbdkms, = -(h • ^fem). 
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We then search for a tensor 

(^ke = (Tle + (Tk£ + (^ke + (^k£ (3-17) 

that must lead to Eq. (3.16) by applying the general relation (3.14). The first 
three terms are simply 

= WeSke - (9fcm • a^m), a^^ = WaSki, (7^ = ''^h^M, (3.18) 

where tUg, Wa and Wh are the energy densities defined in Eqs. (2.19) and (2.23). 
The construction of the magnetostatic contribution is slightly more involved but 
trial and error leads to 

= huh - \h''5M, (3.19) 

where the magnetic induction b = h + m is used mostly as a notational 
abbreviation. As usual, it is understood that the magnetic field in Eq. (3.19) 
is expressed in terms of the magnetization through the magnetostatic equations 
(2.14). A repeated application of these equations establishes the desired relation 

a,aS = -(h-a,m) = Tr. (3.20) 

A notable feature of the derived stress tensor (Jki is that all but the 
magnetostatic contributions are symmetric under exchange of the indices k and 
I. The asymmetry of the last term anticipates the physical fact that the orbital 
angular momentum and the total magnetization are not separately conserved in 
the presence of the magnetostatic interaction, as we shall see shortly. A further 
interesting property is that the preceding construction applies whether or not 
boundaries are present, taking into account that the magnetic induction b is equal 
to the magnetic field h outside the ferromagnetic material where the stress tensor 
reduces to ^ 

(^u = (^k£ = hkhe — -h^Ske, (3.21) 
which is symmetric and satisfies the continuity equation 

dm = (3.22) 

by virtue ofVxh = = V-h. Equation (3.22) is consistent with a vanishing 
Thiele force density outside the material. 



IV. Conservation laws 

The occurrence of ambiguities in the canonical definition of conservation laws 
has already received considerable attention. Slonczweski [16] was apparently the 
first to recognize that the usual definition of linear momentum fails for magnetic 
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bubbles with a nonvanishing winding number. Haldane [17] and Volovik [18] also 
addressed the question from different perspectives. But we believe that a simple 
as well as complete resolution of this issue was given only in some recent work 
[4] where the linear and angular momentum were expressed as moments of the 
topological vorticity (3.1). Since the available studies address strictly 2D and 3D 
models, our aim here is to establish unambiguous conservation laws in the context 
of the quasi-2D geometry appropriate for the description of ferromagnetic films. 

We consider the geometry of Fig. 1 where a film of constant thickness d 
extends to infinity in the {xi,X2) plane and the easy axis is perpendicular to 
the film. Therefore the relevant symmetries are (i) translations in the {xi,X2) 
plane, and (ii) azimuthal rotations around the third axis. We shall show that the 
corresponding conservation laws are the moments 



= y x^,-f3dV, 11 = 1 or 2, (4.1) 

which are related to the linear momentum, and the third component of the total 
angular momentum 

J = e + ^; ^=\j P^lsdV, = Jims- l)dV, (4.2) 

where = + and hence the "orbital" angular momentum £ is also expressed 
as a moment of the vorticity; /x is the total magnetic moment along the easy axis, 
except that we have subtracted the trivial contribution from the ground state so 
that /U is finite and negative. It should be noted that all volume integrals in Eqs. 
(4.1) and (4.2) extend over region I of Fig. 1. 

Although the conservation laws quoted above have a similar appearance with 
those derived for strictly 2D models [4], a proof of their validity is not obvious 
because of potential boundary effects. We consider first the time evolution of the 
moments (4.1), 

4 = y XfiisdV = S3jk J x^djdiakidV, (4.3) 

where we have used the fundamental relation (3.15) applied for i = 3. Notation 
is organized by asserting that Greek indices zv, . . . assume only the two distinct 
values 1 and 2, corresponding to the two spatial coordinates xi and X2, while Latin 
indices i,j,... assume all three values, as usual. We further introduce the 2D 
antisymmetric tensor £^i^, whose elements are en = = £22 and £12 = 1 = —£21, 
and invoke the summation convention for repeated indices without exception. 
Then 

ii^ = e^x J x^d^dtaxidV = e^x J [d^ix^dtaxi) - 5^^dt(Jxi]dV, (4.4) 

where both terms in the integrand are in the form of a total divergence. Since the 
film extends to infinity in the x\ and X2 directions, the first integral vanishes for 
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either = 1 or 2 provided that the magnetization exhibits a reasonable behavior 

at large Xi or X2- By reasonable we mean that the magnetization approaches its 
ground-state value sufficiently fast so that the energy of the configuration is finite. 
Then we may write 

if, = -S/iA / diaxtdV, (4.5) 



where the integrand is also a total divergence but the integral need not vanish 
because the Latin index £ is summed over all three values, £ = 1,2 and 3, and may 
lead to a nonvanishing contribution from the film boundaries, namely 



If, = -^mA 



j axsdxidx2 - J (T\zdxidx2 



xz=d/2 xz ——d/2 



(4.6) 



where the tensor elements uas are evaluated right inside the boundaries and are 
certainly not equal to zero. 

However we may now return to the explicit form of the stress tensor given in 
Section III and apply it for ctas with A 3 to obtain 

= -{d\m ■ dzra) + hxhz, (4.7) 

which must be evaluated at the boundaries of the film where d^m. = on account 
of the unpinned boundary condition (2.22). Hence 

<y\z{x\,X2, xz = ±d/2) = hxbs, (4.8) 

where we further note that the combination of fields hxb3, with A = 1 or 2, is 
continuous across the boundaries thanks to the familiar boundary conditions of 

magnetostatics. The double integrals in Eq. (4.6) may thus be evaluated right 
outside the boundaries where the stress tensor satisfies the continuity equation 
(3.22). An application of the divergence theorem in region II yields 

= y" decrxedV = - J ax3dxidx2, (4.9) 

II X3=d/2 

and a similar relation for region III. The net conclusion is that both integrals in 
Eq. (4.6) vanish and 

4 = 0' (4-10) 

which is the desired result. We shall defer discussion of the interesting physical 
consequences of Eq. (4.10) until a corresponding result is obtained for the angular 
momentum. 

The time evolution of the orbital angular momentum is governed again by the 
fundamental relation (3.15). The analog of Eq. (4.4) now reads 

^ = ^^^A J p'^d^diaxidV = ^e^x J [d^ip'^diaxi) - 2x^d£ax£\dV, (4.11) 
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where the first integral in the last step of Eq. (4.11) vanishes for both i/ = 1 and 
2: 

^ = -ffi/A J x,,de(TxedV = e^^x J Wxi^ - di{x,,axe)]dV. (4.12) 
Recalling that the volume integration extends over region I we write 

J di{xj,axi)dV = J x^axzdxidx2 - J x^ax^dxidx2, (4-13) 

I xz=d/2 xz=-d/2 

where the double integrals may be calculated either above or below the film surfaces 
because the tensor elements axz given by Eq. (4.8) are continuous across the 
boundaries. An argument similar to that used in Eq. (4.9) then leads to 

J x^ax3dxidx2 = - J crxudV, (4-14) 

xz=dl1 II 

and 



x^ax2,dx\dx2 = y oxyti^- (4-15) 

xz = -dl2 III 

Therefore Eq. (4.13) may be rewritten as 

y dt{x,axi)dV = - j axudV - j uxudV, (4.16) 
I II III 

where the right-hand side is symmetric under exchange of the indices v and A 
because the stress tensor is symmetric outside the film. Hence inserting Eq. (4.16) 
in Eq. (4.12) yields a vanishing contribution, because of the contraction with the 
antisymmetric tensor Syx-, and 

(- = j e^xcrxudV. (4.17) 
I 

To summarize, if the magnetostatic interaction were absent, the stress tensor 
would be symmetric in all regions and Eq. (4.17) would lead to a conserved orbital 
angular momentum {£ = 0). In general, using the complete stress tensor given in 
Eqs. (3.18) and (3.19), 

^uxcrxiy = Suxhxbu = Eyxhxmy, (4.18) 

and 

k = j {mih2 - m2hi)dV, (4.19) 
so that the orbital angular momentum is not by itself conserved. 
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Nevertheless a conservation law is obtained by including the total magnetic 
moment fi of Eq. (4.2) whose time derivative is computed by applying directly the 
Landau-Lifshitz equation (2.15) to write 

fi = J m^dV = - y"(m X f)^dV. (4.20) 

Now taking into account the explicit expression for the effective field f of Eq. (2.13) 
we find that the contributions from the anisotropy and bias fields drop out from 
Eq. (4.20) and 

fi = - J {mx Am + mxh)3dV. (4.21) 
To compute the exchange contribution we note that 

J{mx Am)dV = J di{m x d,m)dV = ^"(m x diYn)dSi, (4.22) 

where the surface integral vanishes because of the unpinned boundary condition 
(2.22). Therefore 

ji= - J{m X h)3dV = - J{mih2 - m2hx)dV. (4.23) 

Comparing this result with Eq. (4.19) establishes that 

j = 0, (4.24) 

or that the total angular momentum J = £ + |U is conserved. 

Having thus demonstrated the validity of the conservation laws (4.1) and 
(4.2) we now turn to the discussion of their physical content. We first note 
that these conservation laws are free of all ambiguities even for configurations 
with a nontrivial topological structure. Suffice it to say that the potential 
nondifferentiability of the canonical variables 11 and $ does not affect Eqs. (4.1) 
and (4.2) because they are both expressed in terms of the vorticity which can 
be calculated directly from the magnetization through Eq. (3.1). A detailed 
discussion of this issue may be found in our earlier work within a strictly 2D 
context [4] and applies here with minor modifications. Hence we will simply list 
the important points adapted to the present quasi-2D situation. 

The conserved moments (4.1) are related to the linear momentum p = (^1,^2) 

by 

= ^ixulu-, {Pm' ™} = -^m"^' (4-25) 

where the Poisson-bracket relation establishes that p is indeed the generator 
of translations in the {xi,X2) plane. However p cannot be interpreted as 
ordinary momentum for two related reasons. First, the Poisson bracket of its 
two components, 

{puP2} = -^ndQ, (4.26) 
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does not vanish except for a vanishing winding number. Second, under translations 
in the plane, xi ^ xi + ci and X2 X2 + C2, the moments transform according to 

I^-^ I^ + 47idQc^, (4.27) 

which is a consequence of definition (4.1) and Eq. (3.10). The nontrivial 
transformation of the linear momentum (4.25) implied by Eq. (4.27) is surely an 
unusual property because one would expect the momentum to remain unchanged 
under a rigid translation. Nevertheless the above properties suggest a formal 
analogy with the familiar electron motion in a uniform magnetic field, the role of 
the latter being played here by the winding number. 

Therefore, when Q 7^ 0, a more useful interpretation of the conserved moments 
is obtained through the guiding center coordinates 

which are conserved and transform as (i?i, R2) {Ri + Ci, R2 + C2) under a rigid 
translation in the plane {xi, X2) {xi +ci,X2 + C2). The latter property suggests 
that the 2D vector R = (i?i,i?2) may be interpreted as the mean position of a 
magnetic bubble with Q ^ in a ferromagnetic film, and its conservation implies 
that such a bubble cannot be found in a free translational motion. In other words, 
Q bubbles are always spontaneously pinned or frozen within the ferromagnetic 
medium provided that external perturbations are absent; in analogy with electrons 
undergoing a 2D cyclotron motion in a uniform magnetic film, in the absence of 
electric fields. 

The physical meaning of the orbital angular momentum £ defined in Eq. (4.2) 
is also unusual, for it actually provides a measure of the size of a configuration 
with Q ^ 0. More precisely, one may define a mean squared radius from 

f[ix^-Rif+ix2-R2f]l3dV 
= ^ . = —3- - R^ (4.29) 



/ 



which is directly proportional to £ when the latter is defined with respect to the 
guiding center (R = 0). Note that we use the abbreviated 2D notation R = 
(i?i, R2) and R^ = + i?2- The radius r of Eq. (4.29) plays an important role 
in our theoretical development but does not, in general, coincide with the naive 
radius at which the third component of the magnetization vanishes (ms = 0). 
One should also note that r would be a conserved quantity in the absence of the 
magnetostatic interaction because the orbital angular momentum would then be 
by itself conserved. 

In order to pursue further a meaningful discussion of dynamics, one must 
first ascertain the existence of interesting static solutions of the Landau-Lifshitz 
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equation such as magnetic bubbles with a nonvanishing winding number; an issue 
addressed in the following two sections. We shall return to a more detailed study 
of the implications of the derived conservation laws for dynamics in Section VII. 



V. Virial theorems 

A simple scaling argument due to Derrick [7] leads to a virial relation 
that must be satisfied by any finite-energy static solution of a nonlinear field 
theory. Since Derrick's relation is mainly used in the literature to establish the 
nonexistence of nontrivial static solutions, it is of some interest to demonstrate 
how the present theory evades its potential consequences and leads to the observed 
wealth of magnetic bubbles with practically any winding number [2]. However a 
generalization of the original scaling argument to the present case is not completely 
straightforward, because of the film boundaries, and is given below. 

Static solutions are stationary points of the energy functional W = W{m) 
provided that the constraint = 1 is taken into account. For instance, one may 
use the spherical variables (2.35) to write 

which are the static versions of the Hamilton equations (2.36). In this section, 
we shall neither write out nor solve the above equations explicitly but merely use 
them to derive some general relations. 

For the moment, let us ignore the film boundaries and assume that the 
medium extends to infinity in all directions. We may then apply Derrick's scaling 
argument in a straightforward fashion. Suppose that H = H(x) and $ = $(x) is a 
solution of Eqs. (5.1) with (finite) energy W = We + Wa + Wb + Wm- The energy 
of the configuration H(^x) and $(Cx), where C, is some constant, is then given by 



W{0 = ^We + ^{Wa + Wb + Wm). (5.2) 



By our hypothesis C = 1 is a stationary point of W{() and thus W'{(^ = 1) = or 

We + S{Wa + Wb + Wm) = 0, (5.3) 

which is a virial relation that must be satisfied by any static solution with finite 
energy. Since all pieces of the energy are positive definite, one must conclude from 
Eq. (5.3) that nontrivial static solutions with finite energy do not exist in a 3D 
ferromagnetic continuum without boundaries. 

The preceding derivation of virial relation (5.2) is clearly inapplicable in the 
presence of boundaries. We thus seek to obtain the analog of this relation for the 
film geometry of Fig. 1 by a method that was already employed in the simpler 
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context of Ref. [4] and leads to a series of virial theorems, Derrick's relation being 
the simplest example. An alternative form of the Thiele force density is given by 

5W 5W 

and vanishes for static solutions satisfying Eqs. (5.1). Therefore the stress tensor 
satisfies the continuity equation 

de(TM = (5.5) 

within the ferromagnetic medium. Recalling that the stress tensor satisfies the 
continuity equation outside the film even for time-dependent solutions, see Eq. 
(3.22), we conclude that static solutions satisfy Eq. (5.5) everywhere. 

A series of virial relations may now be derived by taking suitable moments of 
Eq. (5.5) and by a systematic application of the divergence theorem. The simplest 
possibility is 

XjdiakidV = 0, (5.6) 



V 



where the integration extends over some volume V that is left unspecified for the 
moment. The divergence theorem then yields 

J aijdV = j) XjaudSi, (5.7) 

V s 

where we have effected a trivial rearrangement of indices and S is the surface 
surrounding the volume V. It is understood that the region of integration is 
such that the surface S does not cross the film boundaries because of potential 
discontinuities that may render the divergence theorem invalid. 

Thus we proceed with an application of Eq. (5.7) in several steps. First we 
consider the subset of relations obtained by restricting the indices i and j to the 
values 1 or 2. Using our standard convention we write 

I a,^dV = f .^a,,dS„ (6.8) 

V S 

where |U, z/ = 1 or 2, and subsequently apply this relation to each region I, II or 
III separately: 



J a^^dV = -S+^, 



(5.9) 



II 



J a^j,dV = S^^, 



III 
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where 

S^^ = j x^ai^3dxidx2 = J x,,h^jhzdxidx2- (5.10) 

X3=±d/2 xz=±d/2 

Here we have recalled the boundary values of the tensor elements cr^a from Eq. 
(4.8) which are continuous across each boundary for = 1 or 2. In fact, the last 
two equations in (5.9) coincide with Eqs. (4.14) and (4.15) obtained in our earlier 
discussion of conservation laws because the stress tensor satisfies the continuity 
equation outside the film even for time-dependent fields. However the first equation 
in (5.9) applies only to static solutions. An immediate consequence of all three 
equations is the set of relations 

j a^^dV = ^; f^,iy = l or 2, (5.11) 

all volume 

where explicit surface contributions are no longer present. A special case that 
emphasizes the role of the magnetostatic interaction is obtained by contracting 
both sides of Eq. (5.11) with the 2D antisymmetric tensor, 

J ^fn^TnudV = J (mi/i2 - m2hi)dV = 0, (5.12) 

a relation that is consistent with Eqs. (4.19) and (4.23) since both the orbital 
angular momentum £ and the total magnetic moment are time independent in 
a static solution. 

The absence of explicit surface terms in Eq. (5.11) is not surprising because 
scaling arguments of the Derrick variety continue to apply in the xi and X2 
directions. Specifically Eq. (5.11) may be arrived at also by performing the linear 
transformation Xi Cii-"^! +Ci23^2 and X2 — ^ C2i3^i + C223^2 in a static solution and 
by demanding that the resulting energy W = W{Q be stationary at Cn = 1 = C22 
and C12 = = ^21- However the situation is different when one or both indices i,j 
in Eq. (5.7) are equal to 3. 

Actually some useful information on the latter case can be obtained directly 
from the continuity equation (5.5) which is written as 

d„aiu + dsais = (5.13) 

and implies that the double integrals J aizdx\dx2i with i = 1, 2 or 3, are 
independent of xs but may assume different values in regions I, II, or III. In 
fact, all integrals vanish outside the film because they can be calculated at large 
l^sl where the tensor elements vanish. For i — /i — 1 or 2 the integrals vanish also 
inside the film thanks to Eq. (4.9): 

J a^3dxidx2 = 0, // = 1 or 2, (5-14) 
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for any x^. On the other hand, 
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0, l^sl > d/2, 

azzdxidx2 = { (5.15) 
s, -d/2 < X3< d/2, 



where s is constant throughout the film but need not vanish. Collecting the above 
information we may also write 
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assdV = sd, (5.16) 

all volume 

which should be contrasted with Eq. (5.11) where the right-hand side vanishes. 

We have thus derived a number of virial relations that must be satisfied by 
any static solution. We have also gathered sufficient information to make contact 
with relation (5.3) obtained for a strictly 3D medium. Indeed Eqs. (5.11) and 
(5.16) may be combined to yield 
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tTadV = sd, (5.17) 



where the integration extends over all volume, tru = cth + (T22 + cr33 is the trace of 
the stress tensor, s the constant defined from Eq. (5.15), and d the film thickness. 
The trace is calculated by using the explicit expression of the stress tensor from 
Eqs. (3.18) and (3.19): 

3 

trcT = We + 3{wa + Wb) + (h • b) - -b^ (5.18) 

where We,Wa and wi are the exchange, anisotropy and bias energy densities. We 
may further insert in Eq. (5.18) the magnetic induction b = h + m to write 

trcr = We + 3{wa + Wb) - ;Jh^ + 4 ( -^h • m ) - |m^, (5.19) 



.2 J 2 

where the magnetostatic energy density appears both in the form entering Eq. 
(2.25) and that of Eq. (2.32), while in the last term we must set = 1 within 
the film and zero outside. Therefore a more explicit form of Eq. (5.17) reads 

We + 3[Wa + Wb+ (Wm - Wj^^)] = sd, (5.20) 

where we recognize the various pieces of the energy, as in relation (5.3), and Wm^ 
originates in the last term of Eq. (5.19) and is equal to the magnetostatic energy 
of the ground state configuration mo = (0, 0, 1). 

Virial relation (5.20) differs from (5.3) in two significant ways. First, a surface 
term appears in the right-hand side which is entirely due to the film geometry and 
is generally different from zero. Second, the (infinite) magnetostatic energy of 
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the ground state, Wm , is subtracted out. Now implicit in the derivation of (5.3) 
was the assumption that the magnetostatic field vanishes at large distances, in all 
directions, so that the energy is finite. This assumption is clearly false in a 
ferromagnetic film because h = —mo at large xi and X2 and thus both Wm and 
are infinite. Nevertheless the difference — appearing in Eq. (5.20) 
is expected to be finite for reasonable solutions. Furthermore this difference is no 
longer positive definite and is, in fact, negative in the case of magnetic bubbles. 
Indeed the magnetostatic field favors expansion of a domain with magnetization 
opposite to that of the ground state, which is balanced by the exchange, anisotropy 
and bias fields to produce a stable bubble of definite radius [19]. Therefore virial 
relation (5.20), unlike (5.3), does not a priori exclude nontrivial static solutions in 
a ferromagnetic film, irrespectively of the sign of the surface contribution in the 
right-hand side. An explicit example is worked out in the following section where 
both Wm — Wrn ^ and s are negative but Eq. (5.20) is verified. 



VI. The fundamental magnetic bubble 

The construction of static solutions with a nonvanishing winding number is 
an issue of significant practical interest and occupied most of the early studies of 
magnetic bubbles [2, 3] . Because of the long-range nature of the mangetostatic field 
and the related effects of finite film thickness, writing out the static equations (5.1) 
explicitly leads to a rather complex system that is not particularly illuminating. 
Hence the question was addressed through approximate solutions in the limit of a 
large quality factor k [19], variational methods [20], and numerical simulations in 
the important special case of the fundamental [Q = 1) bubble [21]. However, in 
order to illustrate some basic aspects of our theoretical development, we shall need 
some detailed information on the profile of a bubble that is not easily accessible 
from the early work. We have thus decided to recalculate the Q = 1 bubble by a 
numerical method with a simple physical origin. 

Suppose that some initial configuration with a given winding number Q 
evolves according to the Landau-Lifshitz equation (2.37) including dissipation. 
After a sufficiently long time interval precession effects are suppressed and the 
configuration eventually relaxes to a static solution of the Landau-Lifshitz equation 
with the same winding number. Since our aim in this section is only to obtain 
static solutions, the process may be accelerated using Eq. (2.37) with a very large 
dissipation constant A. On introducing the rescaled time variable t — t/X, the 
A — >■ oo limit of Eq. (2.37) reads 

^ + m X (m X f) = 0, m2 = l. (6.1) 
or 

In view of Eq. (2.38) the energy decreases when the configuration evolves according 
to either Eq. (2.37) or its fully dissipative limit (6.1). The advantage of the latter 
is that it suppresses transients and leads to equilibrium with reasonable speed. Of 
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course, the calculated static solution is independent of the details of the initial 
configuration provided that the winding number is kept fixed. Thus the initial 
condition may be chosen more or less at convenience and convergence may be 
improved by incorporating any a priori information on the expected static solution. 

Although the principle of the method is very simple, an efficient solution of 
the initial-value problem posed in the preceding paragraph confronts us with a 
nontrivial numerical task. Thus at every step of the time evolution one must solve 
the Poisson equation (2.27) in order to determine the magnetostatic field h and 
subsequently the effective field f from Eq. (2.13). Calculation of the latter near the 
film boundaries should also take into account the unpinned boundary condition 
(2.22). A detailed description of our numerical algorithm will be given elsewhere 
[22], so the remainder of this section will be devoted to a discussion of the results 
of an explicit calculation of the fundamental magnetic bubble. 

A substantial simplification occurs in the case of the fundamental bubble 
because of its strict axial symmetry; that is, invariance under a simultaneous 
rotation in the {xi,X2) plane and a corresponding azimuthal rotation of the 
magnetization. It is then convenient to use cylindrical coordinates defined from 

xi=pcos(f), X2=psm(f), X3 = z. (6-2) 

A strictly axially symmetric configuration is of the general form 

mi + im2 = (rup + im<^)e*'^, ms = (6-3) 

where the radial {nip), azimuthal (m^) and longitudinal (m^) components are 
functions of only p and z, 

mp = mp{p,z), mct> = m^{p,z), mz = mz{p,z), (6.4) 

while they continue to satisfy the constraint 

m^p + m^ + ml = 1. (6.5) 

The dissipative equation (6.1) becomes effectivelly two-dimensional and a 
significant simplification of the numerical problem results. 

Specifically, when ansatz (6.3) is inserted in Eq. (6.1), the resulting equation 
retains the same form except that the three-component vector m = (mi,m2,m3) 
is formally replaced by {nip, m^, m^) and the effective field f by {fp, f^,fz) with 

fp = Amp «™p + hp, 

= Am^ +hb + hg. 
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where the Laplace operator is reduced to 



92 Id 92 

+ 7 777 + (6-7) 



dp^ p dp dz^ 

hb is the bias field, and hp, and hz are the polar components of the 
magnetostatic field. Actually, the azimuthal component vanishes because 

drup rup dm^ 
V • m = + — + (6.8) 
dp p oz 

and hence the magnetostatic potential is a function of only p and z; = i/^ip, z). 
Therefore 

and the polar components of the magnetic induction are 

bp = hp + mp, b^ = m^, bz = hz + mz. (6.10) 

For future reference we also quote some discrete symmetries of the reduced 
system of equations. First, given a static solution of the form (6.4), the 
configuration 

rripip.z), -m^{p,z), m^ip^z) (6.11) 
is also a solution. Second, the parity relations 

nip (p, z) = -nip (p, -z) , (p, z) = (p, -z) , (p, z) = (p, -z) , 

hp (p, z) = -hp (p, -z) , h^ = Q, hg (p, z) = hg (p, -z) , 

(6.12) 

are compatible with the evolution equation (6.1). In other words, if Eq. (6.1) 
is solved with an initial condition satisfying relations (6.12), the resulting static 
solution will satisfy the same relations. 

To complete the description of strictly axially symmetric configurations we 
return briefly to the conservation laws (4.1) and (4.2). The relevant third 
component of the vorticity reduces to 

1 (a Ti\ 

^' = -p^- ^'-^'^ 

Therefore the winding number calculated from Eq. (3.10) is given by 

1 dm 1 
Q=2j^ -j^dp=-[m,{<^,z)-m,{Q,z)\ = l, (6.14) 

provided that the magnetization approaches its ground-state value rnz = 1 at 
infinity and the value = — 1 at the origin. We further note the trivial fact 
that the moments vanish and the guiding center coincides with the origin of 
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the coordinate system. Finally the orbital angular momentum is computed from 
Eqs. (4.2) and (6.13), 



£ = - dz I —^27Tp'^dp = - / dz (m^ - l)27rpdp, (6.15) 

2 J-d/2 Jo op J-d/2 Jo 



where we have performed a partial integration taking into account that = 1 at 
infinity. We then recognize in the right-hand side of Eq. (6.15) the total magnetic 
moment /x of Eq. (4.2). Hence £ — —jji and 

J = £ + /i = 0. (6.16) 

As expected, the total angular momentum vanishes for a strictly axially symmetric 
configuration. A related fact is that the radius r calculated from Eq. (4.29) with 
R = and I = —ji satisfies the relation 

H = -27rdr^ (6.17) 

which could also be arrived at by considering a crude model of a bubble where the 
magnetization points toward the north pole, m = (0,0, 1), for p > r and toward 
the south pole, m = (0, 0, —1), for p < r. 

Now, if Eq. (6.1) is solved for an initial condition with strict axial symmetry 
and winding number Q — 1, it will eventually lead to a static solution with the 
same symmetry and winding number. A simple choice of the initial configuration 
is given by the two-parameter family 

nip = 0, m0 = ±sech«, = tanhu, (6.18) 

with 

u = Hp/po) + {p- po)/So, (6.19) 

which coincides with the variational ansatz employed by DeBonte [20] treating the 
constants po and So as variational parameters. The constant po is the naive radius 
of the bubble, i.e., the radius at which the third component of the magnetization 
vanishes, while both po and Sq provide a measure of the wall width Syj in a picture 
where the bubble is viewed as a curved domain wall: 



1 du 
Sw dp 



P=Po 



^ + 1. (6.20) 

Po Oo 



On this occasion we recall that the width of an ideal (straight) domain wall in an 
infinite medium is 

^^ = \ or -\= (6.21) 



K 



in the original or rationalized units, respectively (see Section II). Needless to say, 
for our purposes the constants po and 5o need not be determined variationally 
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because the relaxation algorithm should lead to the true static solution for any- 
choice of these parameters. However convergence may be accelerated when 
configuration (6.18) is as close as possible to the true bubble. 

The description of the initial ansatz is completed noting that the ± freedom 
in Eq. (6.18) reflects the discrete symmetry (6.11). The specific choice of sign in 
will be referred to as the polarity of the bubble, the winding number being the 
same (Q = 1) for either polarity. Finally configuration (6.18) is independent of z 
and trivially satisfies the parity relations (6.12). Therefore the anticipated static 
solution will satisfy the same relations, even though it will develop a nontrivial z 
dependence. 

At this point, one must specify the true parameters of the problem, namely 
the quality factor k, the bias field hi, and the film thickness d. We have aimed 
at providing an illustration where the bubble radius is roughly equal to the film 
thickness and have thus arrived at the specific values (in rationalized units) 

1 6 

K = 2, /lb = 0.32, d=16A„ = ^, (6.22) 

which belong to a parameter regime that is thought to be ideal for the formation of 
magnetic bubbles [19]. A possible choice of the variational parameters in the initial 
ansatz (6.19) is accordingly given by po = 18 and do — l.lA^ but it is certainly 
not unique. We finally mention that in all of the ensuing graphical illustrations of 
the fundamental bubble we invoke a slight departure from the rationalized physical 
units introduced in Section II and used throughout our theoretical development. 
Thus distances will now be measured in units of the ideal domain wall width 
A^ = I/a//^ = 1/a/2 in order to emphasize the wall structure of the calculated 
bubble. For instance, the film thickness will appear as d — 16. 

The calculated fundamental magnetic bubble is illustrated in several ways. 
We mostly describe a Q = 1 bubble with positive polarity, originating in the 
initial ansatz (6.18) with the upper sign in m^, the results for negative polarity 
being inferred from the discrete symmetry (6.11). Thus in Fig. 2 we display 
the dependence of the magnetization on the radial distance p at the film center 
(z = 0) and near the upper boundary {z = d/2); the p dependence near the lower 
boundary {z = —d/2) may be obtained from the parity relations (6.12). The 
corresponding results for the magnetic induction are shown in Fig. 3. One should 
keep in mind that the magnetostatic field extends beyond the film boundaries, but 
the calculated values will not be discussed further in the present paper. 

Some important general features of the fundamental bubble are already 
apparent in Fig. 2. If we view the bubble as a curved domain wall, the wall 
is purely Bloch at the film center {nip = 0) and nearly Nccl at the boundaries 
where the radial component nip achieves significant values while the azimuthal 
component is small. A better view of the situation is obtained by plotting the 
projection of the magnetization vector m on the (xi, X2) plane in Fig. 4; whereas 
Fig. 5 illustrates the projection on a plane that contains the easy axis, which 
is chosen to be the (xijXs) plane without loss of generality thanks to the axial 
symmetry. 
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The case of a Q = 1 bubble with negative polarity may be inferred from Eq. 
(6.11) applied, for example, to Fig. 4. The reflection —rriff, will reverse the 

sense of circulation of the magnetization at the film center but will not significantly 
affect the picture at the boundaries where is small. 

Having thus provided an overall view of the fundamental bubble, we now turn 
to the description of some important details. The simplest way to analyze the fine 
structure of the bubble is by recalling the concept of a vortex line introduced in 
Section III. We first restate the condition of strict axial symmetry in terms of the 
spherical variables (2.33): 

Q = 9{p,z), $ = + x(p,^), (6.23) 

where the functions 9 and x cire independent of the angle (f) and are related to the 
polar components of the magnetization by 

cos 9 — rriz, x = arctan(m<^/mp). (6.24) 

Therefore a vortex line is equivalently defined as the intersection of the two surfaces 

mz{p, z) ^ 1713, + x(p, ^) = 0o, (6.25) 

where ms and 0o ^ire arbitrary constants in the intervals [—1,1] and [0, 27r], 
respectively. 

The first relation in (6.25) defines a curve in the (p, z) plane, illustrated in Fig. 
6 for three typical values of ms, and the surface obtained by a simple revolution of 
the curve around the third axis has the shape of a barrel. Of special interest is the 
case 777-3 = which may be used to define a (naive) radius of the bubble po — po{z) 
as the radius of the circular intersection of the barrel with the {xi,X2) plane at 
altitude z. We use the same symbol for the naive radius as for the variational 
parameter po in Eq. (6.19) because the two coincide within the initial ansatz. The 
current definition of the radius is especially useful at the film boundaries where po 
is the distance from the center of the bubble at which a sharp change in contrast 
takes place (see Fig. 4) that may be detected experimentally. In our numerical 
example we found that pq{z = ±d/2) = 15. 9 A^^, which should be compared with 
the value at the film center po{z = 0) = 16. 3 A^,, thus providing a measure of 
bubble bulging [19]. We also quote the average value of the naive radius 

1 rd/2 

Po = / Po{z)dz = 16.15A„ (6.26) 

" J-d/2 

and compare it with the radius r defined by Eq. (4.29) and related to the total 
magnetic moment by Eq. (6.17): 

r = 16.21A^. (6.27) 

The radius r appears naturally within the theoretical development, as will become 
evident in the discussion of skew defiection in Section VII, whereas the naive radius 
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Po is closer to what is actually measured in an experiment. Therefore the observed 
proximity of the numerical values quoted in Eqs. (6.26) and (6.27) is of special 
significance. Note incidentally that the calculated bubble radius for the specific 
parameters (6.22) is approximately equal to the film thickness. 

Next we consider the second relation in (6.25). We will not attempt to draw 
the corresponding surfaces, for various values of the constant (f)o, but examine 
directly their intersections (vortex lines) with a barrel at given ms- Thus in Fig. 
6 we also display the ^-dependence of x along a vortex line; that is, we plot 
the function x{p{z; ms) , z) = x(2;;m3) where p{z\m^) is the root of the algebraic 
equation 'mz{p, z) = which depends on z and the particular value of m^. At 
777,3 = 0, the root p^z^ms — 0) reduces to the naive radius po{z) discussed in the 
preceding paragraph. Actually Fig. 6 shows the ^-dependence of x along a vortex 
line only for ma = 0, but our numerical simulation furnished values for x that are 
virtually indistinguishable when ms varies in the region \mz\ < 1/2. Nevertheless 
departures from such a universal behavior occur for \mz\ > 1/2. 

The preceding findings are best summarized by the sketch of a typical vortex 
line given in Fig. 7. Since the sum (f) + x(p(2;m3),2) must remain equal to a 
constant cj)Q, a vortex line originating at a point A of the upper boundary, i.e., at 
given 00 and ma, proceeds downward along the surface of the barrel at the same 
time twisting by an amount determined by the variation of x illustrated in the 
lower part of Fig. 6. The vortex line eventually terminates at the lower boundary, 
at a point B, having suffered a total twist = —Ax = x(z = d/2; nis) — xi^ = 
—d/2; ms) which generally depends on m3. For the specific example m3 = shown 
in Fig. 6 the calculated total twist is = 158°. This value remains practically 
the same in the range jmsj < 1/2 but deviations do occur for |m3| > 1/2. Finally 
we note that the vortex lines twist around the surface of a barrel counter-clockwise 
for the Q = 1 bubble with positive polarity considered in our illustrations; the twist 
takes place clockwise for the Q — 1 bubble with negative polarity obtained through 
the discrete symmetry (6.11). 

Returning to the topological classification discussed in Section III, we note 
that all vortex lines in the fundamental magnetic bubble terminate at the film 
boundaries and do not tangle with each other. Hence a definition of a Hopf index 
is not possible, as expected for a film of finite thickness. On the other hand, the 
flux of vorticity through the plane {xi, X2) is equal to Att for all z, thus leading to 
a unit winding number. 

We also return to the virial relations derived in Section V and comment on 
the manner they are satisfied in our explicit calculation of the fundamental bubble. 
Because of the strict axial symmetry the surface integrals S^^^ of Eq. (5.10) reduce 
to 




(6.28) 



where 6,,i, is the 2D Kronecker delta and 




^ Jo 

X3 = ±d/2 



(6.29) 
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Hence relations (5.9) applied ioi fj, ^ v simplify to 



J ai2dV = = y" (^2idV, (6.30) 

which are valid in each region I, II or III separately. In fact, using the explicit 
expression of the stress tensor for an axially symmetric configuration, Eqs. (6.30) 
reduce to the single equation 

j hpm^dV = 0, (6.31) 

which is automatically satisfied thanks to the parity relations (6.12). Furthermore 
Eqs. (5.9) applied for = = 1 or 2 yield for the various regions 

j aiidV = J (T22dV = S+-S-, 
I I 

J a^^dV = J (j22dV = 



(6.32) 



II II 



j (Tiidy = j a22dV = S-, 
III III 



as well as 



J GiidV = = y (y22dV, (6.33) 



when the integration extends over all volume. The above relations were verified 
explicitly in our numerical simulation and were thus used to test its validity. We 
further verified the Derrick-like relation (5.20) in our specific numerical example 

where We = 4134, Wa = 6008, Wi, = 11952, Wm - = -22573, and sd = 

—9713; here energy is measured in the rationalized units introduced in Section II. 



VII. Skew deflection 

We are thus ready to study the main dynamical question posed in the 
Introduction. A static bubble with winding number Q initially located at, say, 
the origin of the coordinate system is subjected to an external magnetic field 

hext = (0,0, /text), /text = ^ext(x,t), (7.1) 

that points along the easy axis and its strength is some prescribed function of 
position and time. Our task is to determine the response of the bubble to such 
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an external probe. In the absence of dissipation (A = 0) the relevant dynamical 
equation is Eq. (2.15) extended according to 



f^f+hext, (7.2) 

in order to include the effect of the applied field (7.1) which is turned on at t = 0. 
One must then solve the resulting equation with an initial condition provided by 
the static bubble and calculate the magnetization m = m(x, t) at all later times. 

However a great deal can be learned without actually solving this initial- 
value problem thanks to the special nature of the conservation laws derived in 
Section IV. Since the position of the guiding center R = (i?i,i?2) is conserved 
in the absence of dissipation and external fields other than a uniform bias field, 
examining the rate at which R changes in the presence of the field (7.1) should 
yield direct information on the response of the bubble. Hence the vorticity now 
obeys the relation 

ii = Sijk [djdeake - dj (hext • ^fem)] , (7.3) 

which is Eq. (3.11) or (3.15) extended according to Eq. (7.2). In particular, the 
evolution of the third component of the vorticity is governed by 

73 = £u\ [di^d^axe - d„ (hext • ^Am] , (7.4) 

where we have returned to the 2D notation for Greek indices, as in Section IV, 
except for the Latin index £ that is summed over all three values. The evolution 
of the moments 1^ of Eq. (4.1) is then given by 

-^M = y" £i^xXij.[d^di(Tx£ - 5i.(hext • dxm]dV, (7.5) 

where the contribution of the first term may be shown to vanish by reasoning 
completely analogous to that used in the derivation of the conservation laws in 
Section IV. Implicit in the above statement is the assumption that the applied 
field (7.1) does not affect significantly the configuration of the bubble at large 
distances or, equivalently, it does not affect the ground state of the ferromagnet. 
We shall return to this assumption later in this section. Thus the evolution of the 
moments in the presence of the applied field is governed by 

ii^ = - J £uXX^d^{he^tdxm3)dV = J e^^he^tdumsdV, (7.6) 

where we have taken into account that the field (7.1) points in the third direction 
and have also performed an elementary partial integration. On the assumption 
that 777-3 1 sufficiently fast at spatial infinity, one may perform a further partial 
integration to write 

= {ef^Me^t){ms - l)dV. (7.7) 
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Also recall that the winding number is conserved even in the presence of the applied 
field provided that the latter does not destroy the ground state of the ferromagnet. 
Therefore the drift velocity of the bubble may be inferred from Eqs. (4.28) and 
(7.7): 

V^, = R^ = J ie^Me.t)im3 - l)dV. (7.8) 

This result for the drift velocity is not completely explicit because the third 
component of the magnetization appearing under the integral sign must still be 
determined through a detailed solution of the initial-value problem described in 
the introductory paragraphs of this section. However Eq. (7.8) already contains 
the essential information concerning the experimentally observed skew deflection 
of magnetic bubbles, for it suggests that the drift velocity will acquire a significant 
component mainly in a direction perpendicular to the gradient of the applied field. 

In order to appreciate the physical content of Eq. (7.8) the applied field is 
written as 

he^t = gxi, g = g{yi,t), (7.9) 

where the "gradient" g may still be a function of position and time. The two 
components of the drift velocity (7.8) are given equivalently by 

Vi = J {x,d2g){ms - l)dV, V2 = ^ J {g + x^d^g){m^ - l)dV. 

(7.10) 

The field is now restricted to the physically interesting situation where the gradient 
is nearly spatially uniform, i.e. g ~ g{t), over a large region surrounding the 
bubble and drops to zero outside that region. Under such conditions all implicit 
assumptions made in deriving Eqs. (7.10) are satisfied. In particular, all partial 
integrations performed on the assumption that the field does not significantly alter 
the behavior of the bubble at large distances are justified. Nevertheless it is clear 
from Eqs. (7.10) that an especially transparent result would be obtained in the 
ideal limit where the gradient g is spatially uniform everywhere. Then 

» = ^'=»' ''^ = 4^' ("1) 

where ^ is the total magnetic moment defined in Eq. (4.2). 

A completely uniform gradient would imply an infinite field at xi — > —00 
opposing the magnetization in its ground state m = (0,0,1). Therefore, in the 
presence of some dissipation, the magnetization would align with the applied 
field almost immediately after the field is turned on and assume the value 
m = (0, 0, —1) far in the left plane. Such an instance would destroy the original 
topological structure of the bubble and obscure the question of skew defiection. 
Nevertheless the preceding criticism does not apply in the case of vanishing 
dissipation considered so far because the magnetization would then precess wildly 
around the easy axis far in the left plane but never align with the external field. 
Hence the uniform limit considered in Eq. (7.11) is mathematically meaningful in 
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the absence of dissipation and provides the clearest, albeit idealized, illustration 
of the main theme discussed in this paper. 

Thus the guiding center of a bubble with Q 7^ moves in a direction 
perpendicular to the applied gradient, in analogy with the familiar Hall motion 
of an electric charge in a uniform magnetic field (the analog of the winding 
number) and a uniform electric field (the analog of the uniform magnetic-field 
gradient). Furthermore the drift velocity (7.11) is expressed in terms of quantities 
with a simple physical meaning. However this expression for the drift velocity 
is not completely explicit because the total moment is not conserved during the 
application of the gradient and thus acquires some time dependence fj, = n{t) that 
can be determined only through a detailed solution of the initial-value problem. 
The moment would be conserved if the magnetostatic field were absent. Indeed 
the orbital angular momentum £ and the total magnetic moment fi would then 
be separately conserved in the absence of the applied field and, while the latter 
violates conservation of £ because it breaks rotational symmetry, it does not afi^ect 
H because it points in the third direction. Under such conditions Eq. (7.11) 
provides an explicit expression for the drift velocity since the conserved moment 
// could then be calculated from the initial configuration of the (static) bubble. 
This situation occurs in the case of the 2D isotropic Heisenberg model where an 
analytical result for the drift velocity was given in Ref. [4] and was later verified 
by a numerical simulation in Ref. [5]. 

Returning to the realistic case where the magnetostatic field is not negligible, 
we note that the total moment fi may still be calculated from the profile of the 
static bubble during the initial stages of the process. If we further restrict our 
attention to the fundamental bubble calculated in Section VI, the moment is 
related to the bubble radius through Eq. (6.17) and the initial drift velocity 
may be written as 

g = l; V,^0, V2 = -\gr^ (7.12) 

where we may substitute the numerical value for the bubble radius given in Eq. 
(6.27) which is approximately equal to the naive radius measured in an experiment. 
Now applying Eq. (7.12) for the speed V = IV2I yields 

which is golden rule (1.1) with a defiection angle 5 = 90°, as is appropriate in the 
absence of dissipation, and a winding number Q = 1. 

However one should keep in mind that the above result provides only a 
partial verification of the golden rule for two reasons. First, Eq. (7.13) is in 
general violated during the late stages of the process because neither /x nor r 
are conserved; in particular, the relation // = — 27r(ir^ is strictly valid only for a 
static bubble. Second, one must examine the extent to which (1.1) is valid in 
the presence of dissipation when the defiection angle is no longer equal to 90°. 
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The last statement follows from the simple physical fact that dissipation induces a 
tendency for alignment of the magnetization with the external field, which drives 
the bubble also toward the left half plane where the field points along the negative 
third direction. 

In order to study the effect of dissipation more precisely we return to Eq. 
(2.37) and extend it according to Eq. (7.2) to write 

m + (mxG)=0, = 1, (7.14) 

where the effective field G is given by 

G = Ai(f + hext) + A2[m X (f + hext)], 

- TTa^' - TTa^- 

Since Eq. (7.14) is formally identical to the dissipationless equation, with the 
replacement f — > G, the time derivative of the vorticity may be inferred from Eq. 
(3.11) with the same replacement: 

ii = -SijkdjiG ■ dkm). (7.16) 

Substitution of the field G from Eq. (7.15) leads to 

ii = eijkdj{Xidiake - Ai(hext ■ dkUi) - A2[m x (f + hext)] • dkm}, (7.17) 

a result that may be used to study the time evolution of the guiding center in a 
manner analogous to our earlier discussion in the absence of dissipation. The drift 
velocity is now given by 

V,, = £^ J [Aihext + A2[m X (f + hext)] • d.mdV, (7.18) 

which reduces to Eq. (7.8) at vanishing dissipation. 

The above result is highly implicit in that the magnetization in the right- 
hand side must still be determined through a detailed solution of the initial value 
problem. However some explicit information can be extracted from Eq. (7.18) for 
the early stages of the bubble motion. The magnetization may then be calculated 
from the static profile of the bubble for which m x f = 0. Therefore the initial 
drift velocity is given by 

= £^ J [Aihext + A2(m X hext)] " d^mdV, (7.19) 

where it is understood that the magnetization is that of a static bubble with 
winding number Q. Taking into account that the field points in the third direction 
and performing a partial integration in the first term yields the equivalent relation 

= -^^^ J [Ai(9^/iext)(w3 - 1) + X2he^t{mid^m2 - m2d^mi)]dV. (7.20) 
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At this point, one should recaU the assumptions on the gradient g = ^(x, t) 
discussed following Eq. (7.9) which are especially important in the presence 
of dissipation. Specifically the gradient must vanish outside a large region 
surrounding the bubble, for otherwise both the ground state and the topological 
structure of the bubble would be significantly altered. Yet the specific choice of 
the gradient at large distances will certainly affect the long-time behavior of the 
bubble but should not be crucial during the early motion. Therefore it is still 
meaningful to approximate the initial drift velocity by inserting in Eq. (7.20) the 
applied field /lext — gxi with a gradient g = g{t) that is spatially uniform. If 
we further restrict Eq. (7.20) to the fundamental magnetic bubble calculated in 
Section VI we find that 

where A is the dissipation constant, d is the film thickness, and the constants v 
and c are given by 

jjL^ I {niz - l)dV, v = ^ I {ml + rni)dV, 



2 



(7.22) 



Here fi is the total moment, v is essentially the anisotropy energy, and c vanishes 
on account of the parity relations (6.12). Therefore our final result for the initial 
drift velocity in the presence of dissipation is 

V. = -X.£-^, V,=A.£i. (7.23) 

Since the total moment ^ is always negative and the constant u positive, the 
guiding center moves off in the lower left plane with an initial defiection angle 5 
with respect to the negative xi axis given by 

Vo 1 1 
tan 5 = — = — or sm.5 = :, (7.24) 



where the coefficient 



7/ =-- = 0.08 (7.25) 

A* 



is calculated from Eqs. (7.22) using as input the static bubble. The specific 
numerical value quoted above corresponds to the specific choice of parameters 
made in Section VI. Finally we calculate the speed 
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and relate it to the bubble radius through Eq. (6.17) and the deflection angle 
through Eq. (7.24): 

^—sm5=^-^. (7.27) 

At vanishing dissipation (A = 0) the deflection angle (7.24) becomes 5 = 90° and 
relation (7.27) reduces to (7.13). 

Relation (7.27) establishes contact with the semi-empirical golden rule (1.1) 
in the important special case of the fundamental magnetic bubble. For small 
values of the dissipation constant A encountered in practice [2] the right-hand 
side of Eq. (7.27) is well approximated by unity, to within terms of order A^, 
and is thus consistent with Eq. (1.1) applied for Q = 1. However this is again 
only a partial veriflcation of the golden rule because Eq. (7.27) is strictly valid 
only for the initial drift velocity. A complete verification would require first to 
ascertain that the bubble eventually reaches a steady state, namely a state with 
constant velocity and radius. Such a question could be addressed by a direct 
numerical solution of the initial- value problem posed in the first paragraph of this 
section. This numerical task is in several respects similar to the solution of the fully 
dissipative equation (6.1) described in Section VI, except for a technical difference 
that might prove crucial in practice. Because the applied field breaks rotational 
invariance the bubble looses its strict axial symmetry during skew deflection and 
thus leads to a 3D numerical simulation that is beyond our current capabilities. 

Nevertheless the question was addressed and answered within a strictly 2D 
Skyrme model [6] which also leads to Eqs. (7.24) and (7.27) for the initial drift 
velocity. However it was found through an explicit numerical solution of the initial- 
value problem that a sharp transient period exists during which the deflection 
angle departs rapidly and signiflcantly from Eq. (7.24). The transient period is 
followed by an intermediate regime where the deflection angle reaches a more or 
less constant value and the golden rule is verifled in a rough manner. But a true 
steady state is never achieved and the finer predictions of the golden rule are not 
sustained. In particular, the long-time behavior of the bubble is sensitive to the 
details of the gradient at large distances. 

The results from the 2D Skyrme model [6] could be used as a guide for future 
numerical investigations of the realistic quasi-2D model studied in the present 
paper. We thus turn to a summary of our main conclusions given in the following 
section. 



VIII. Concluding remarks 

We believe to have provided a clear illustration of an important link that 
exists between the topological complexity of ferromagnetic structures and their 
dynamics. The most direct manifestation of such a link is the construction 
of unambiguous conservation laws as moments of the topological vorticity. 
The special dynamical features of magnetic bubbles become transparent and 
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are formally related to more familiar situations such as the Hall effect of 
electrodynamics or the Magnus effect of fluid dynamics. 

Our work has also revealed that some of the quantitative predictions of the 
early studies must be interpreted with caution. In particular, the golden rule is 
valid in its gross features but not in its details. Hence there exists room for further 
development of the dynamical theory of magnetic bubbles. Numerical simulations 
along the lines of those performed within the strictly 2D Skyrme model [6] could 
prove feasible and provide important hints concerning the remaining questions. 
There has also been some speculation to the effect that the dynamics might simplify 
for hard {\Q\ 3> 1) bubbles, in analogy with the adiabatic dynamics of electric 
charges in strong magnetic fields [5]. The scnii-cmpirical golden rule might then 
prove to be exact in the extreme large-Q limit and could possibly be corrected 
through a systematic adiabatic perturbation theory at finite Q. 

In this paper we have confined our attention to the response of a bubble to an 
externally applied magnetic-field gradient. However a field gradient is intrinsically 
present also in the problem of two or more interacting magnetic bubbles. Thus two 
interacting bubbles with winding numbers of the same sign are expected to orbit 
around each other, in analogy with the 2D motion of two electrons in a iniiform 
magnetic field or two vortices in an ordinary fiuid. Similarly two bubbles with 
opposite winding numbers (e.g., Q = 1 and Q = —1) should move in formation 
along roughly parallel lines, also in analogy with an electron-positron pair in a 
uniform magnetic field or a vortex-antivortex pair in a fluid. These expectations 
were confirmed through numerical simulations in the Skyrme model [6] and should 
be possible to establish both theoretically and experimentally in real ferromagnetic 
films. Analogous results are expected for interacting Abrikosov vortices in a 
superconductor [23]. 

Finally we return briefiy to the possibility of genuinely 3D magnetic solitons 
alluded to in Section III. The experimentally observed Bloch points are 3D 
topological defects whose dynamics has not yet been studied within the present 
framework. Furthermore theoretical arguments for the existence of magnetic 
vortex rings with a nonvanishing Hopf index [4, 15] have not yet been concluded 
to a definite calculation that would provide the necessary background for a 
corresponding experimental search in the bulk of the ferromagnetic medium. 
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Figure captions 



Figure 1: Conventions concerning a ferromagnetic film of thickness d. The film 
extends to infinity in the xi and X2 directions. 

Figure 2: The calculated magnetization for the fundamental {Q = 1) magnetic 
bubble with positive polarity and parameters specified by Eq. (6.22). Results are 
given at the film center {z = 0) and at the upper boundary {z = d/2), whereas the 
corresponding results at the lower boundary [z — —d/2) may be inferred from the 
parity relations (6.12). Here and in all subsequent graphical illustrations distance 
is measured in units of the ideal wall width A^^ = 1/ v^- 

Figure 3: The calculated magnetic induction for the Q = 1 bubble; see caption 
of Fig. 2 for further explanations. 

Figure 4: Illustration of the Q = 1 bubble with positive polarity through the 
projection of the magnetization vector field m on the {xi,X2) plane. The bubble 
is Bloch-like at the film center (xs = 0) and Neel-like near the boundaries (^3 = 
±d/2). 

Figure 5: Another view of the Q = 1 bubble through the projection of m on the 
(xijXs) plane. 

Figure 6: The calculated curves mg{p,z) = ms for = 0, ±1/2 (upper entry) 
and the 2;-dependence of x along a vortex line (lower entry) for the Q = 1 bubble 
with positive polarity; see the text for further explanations. 

Figure 7: Sketch of a typical vortex line for a Q = 1 bubble with positive polarity. 
The sense of twist is reversed for a Q = 1 bubble with negative polarity. 
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